A Langevin analysis of fundamental noise limits in Coherent Anti-Stokes Raman 

Spectroscopy 
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We use a Langevin approach to analyze the quantum noise in Coherent Anti-Stokes Raman 
Spectroscopy (CARS) in several experimental scenarios: with continuous wave input fields acting 
simultaneously and with fast sequential pulsed lasers where one field scatters ofi' the coherence 
generated by other fields; and for interactions within a cavity and in free space. In all the cases, 
the signal as well as the quantum noise due to spontaneous decay and decoherence in the medium 
are shown to be described by the same general expression. Our theory in particular shows that 
for short interaction times, the medium noise is not important and the efficiency is limited only by 
the intrinsic quantum nature of the photon. We obtain fully analytic results without making an 
adiabatic approximation, the fluctuations of the medium and the fields are self solved consistently. 

PACS numbers: 42.65.Dr, 05.10.Gg, 05.40.Ca, 87.64.Je 



I. INTRODUCTION 

The field of Raman spectroscopy is a very mature one 
that boasts a vast literature which over time has de- 
veloped and explored countless ingenious improvements 
and techniques to tweak out better and stronger signals. 
Much of the motivation for this lies in the fact that Ra- 
man scattering inevitably involves the structure of the 
scattering medium, because the incoming and outgoing 
signals differ by the frequency of some internal mode of 
the constituent molecules or atoms, thereby making it an 
invaluable tool for spectroscopy. 

The weak signals associated with spontaneous Raman 
scattering were long overcome by stimulated scattering 
through non- linear interactions. Signal has been en- 
hanced through resonance with some specific natural 
modes of the molecules. In particular the technique 
of Coherent Anti-Stokes Raman Scattering (CARS) has 
emerged as one of the most useful Raman technique; 
it involves two incoming fields that create coherence in 
the medium which thereby enhances the scattered sig- 
nal. The generated field being anti-Stokes eliminates flu- 
orescence problems, and resonance enhancement is also 
usually possible. 

One way to improve on existing CARS techniques has 
been recently proposed which uses sequential fem- 
tosecond pulses whereby maximal coherence is created 
via adaptive algorithms with one set of lasers and then 
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the medium is probed by another laser. The technique 
called FAST (for Femtosecond Adaptive Spectroscopic 
Techniques) CARS has been applied in preliminary ex- 
periments 13 and further promises to increase the CARS 
Raman signal. Experimental work is underway at sev- 
eral laboratories. The method offers hope for developing 
a way of detecting (in real time) dangerous biological 
spores like anthrax. 

For any spectroscopic method its relevancy and effec- 
tiveness is eventually defined by the signal/noise ratio. 
A brief review through the literature [3, Q shows that 
there have been methods developed to overcome almost 
every possible laboratory source of noise be it be unsta- 
ble lasers or non- resonant background signals, perhaps 
not all with the same technique but the point is that 
such techniques exist. But at the end there is still the in- 
herent quantum noise of the system which is unavoidable. 
There are two sources for this. The first is truly funda- 
mental in the sense that it represents the lower limit of 
signal detection, this is the shot noise which is associated 
with Poisson statistics of lasers used far above threshold. 
The second arises from the spontaneous emission and the 
decoherence associated with excited atoms and molecules 
created during all Raman scattering processes. Often this 
second source of noise, that we will refer to as medium 
noise, is much larger than the shot noise and therefore is 
the limiting quantum noise. 

Thus our goal in this paper is to undertake a compre- 
hensive study using general Langevin methods to analyze 
the quantum noise in CARS and FAST CARS. Our ap- 
proach and considerations have the advantage of being 
general and may easily be used to study other coherent 
Raman processes. We study several experimental scenar- 
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FIG. 1: CARS scheme with three input fields of frequency vi, 
V2 and uz and a generated signal field V4, = ui -\- us — V2. The 
radiative damping rate from level b to level c is Fj, and the 
dephasing rate of the coherence between levels b and c is F. 



ios allowing for interactions within or without a cavity 
and for pulsed or continuous wave input fields. We find 
a remarkable similarity in behavior of all these cases sug- 
gesting that our results probably have a broader validity 
beyond the configurations we consider. In particular, we 
find that the medium (solvent) noise is not important for 
FAST CARS whose efficiency is limited by the intrinsic 
quantum character of the photon. 

In sections II and III we define the problem in terms 
third order non-linear interactions and Langevin equa- 
tions, and in section IV we describe our model in detail, 
laying out the equations and assumptions that we use. 
In Sec. V we consider the experiments using sequential 
femtosecond pulses and in Sec. VI we discuss experi- 
ments involving concurrent continuous wave input fields. 
Then in Sec. VII we present a detailed discussion of our 
results and their implications, and we provide numerical 
estimates based on realistic experimental parameters. 



II. NONLINEAR INTERACTION 



all non-resonant terms. For the case of FAST CARS 
where vi and V2 pulses interact with the atoms before 
1^3, the non- resonant terms are simply not present. How- 
ever when all the fields are present simultaneously the 
permutations of the fields lead to 24 terms in the pertur- 
bation expansion which all contribute at the same order 
; in that case we have in effect neglected the 20 remain- 
ing terms which have non-resonant combinations of the 
fields and the atomic levels. If we are close to resonance 
this is certainly justified, in fact we will assume perfect 
two photon Raman resonance in this paper. 

The Raman coupling constants are defined by R12 ~ 
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{jlfj ■ei){ilji -62) iftfj ■£2)(Mj» -ei) 



{ujjc + V2) 



(2) 



and likewise for i?43. The e are unit polarization of the 
fields, the dipole moments, j labels the intermediate 
states and £1 = \J (hui) / {2€oepV) are the field quanti- 
zation factors over volume in a medium of dielectric 
constant ep taken to be about the same for all laser fields. 
The decoherence rate between levels b and c is denoted 
by r and the depopulation rate from level b to level c by 
r^; damping to other levels is assumed negligible. 

It will help to establish a relation with macroscopic 
quantities. The classical field amplitudes are Ei = EiUi 
for coherent state expectations ai = {cLi). The third- 
order susceptibility is defined to be 



with expected population difference between the two lev- 
els N = Nf, — Nc- The Maxwell equation for the slowly 
varying amplitude of the generated field is 



ot oz J 2eoep 



(4) 



with the velocity in the medium v = cj ^ft^. 



We will consider stimulated Raman scattering involv- 
ing three input fields ai, 02, 03 and one generated field 04 
with corresponding frequencies 1^^=1,2, 3, 4- The interac- 
tion Hamiltonian is written in terms of the two Raman- 
coupled atomic or molecular levels |6) and |c), of energies 
hioi, and hiOci ^^id lowering operator a = |c)(6| 



V = 'h[R 



112O1O2 



H.a. 



(1) 



Rere and elsewhere the time arguments will he suppressed 
where obvious in order to reduce clutter. The exponen- 
tial factor signifies that all operators are in the interac- 
tion picture, its argument A = w^c + 1^21 — ^bc + 1^34 
with LJiic — ~ ^c- These particular combination of 
atomic and field frequencies in A implies that the Hamil- 
tonian contains only Raman resonant terms and excludes 



III. LANGEVIN EQUATIONS OF MOTION 

The three input fields vi , V2 and z/3 are typically strong 
laser fields far above threshold, we will treat them as clas- 
sical fields and replace 01(2,3) ~^ 0^1(2,3)- Considerable 
simplification is achieved by defining 



Ri2aia2 + R43ala4{t) 



Ri2aia2 
9 



g = \Ri3,al 



(5) 



The normalization g has been chosen such that A(t) sat- 
isfies the boson commutation rules. Thus we can treat 
A{t) like a photon field operator. The parameters g and 
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^ being dependent only on the input fields will treated as 
constants in our calculations. The Heisenberg-Langevin 
equations of motion for the atomic operators are obtained 
from the interaction Hamiltionian: 



TbOb - ig 



A^a-a'^A 



+ Fb 
Fc. 



(6) 



The coherence operator here differs from that in Eq. 
by a phase factor a{t)e^^^* ~> o-(t), and ab(c) are the 
population operators for the atomic levels. The F^s are 
the Langevin noise operators. They have vanishing first 
moments {F{t)) — 0, since the entropy of the system 
cannot be lowered by noise The second moments 
are taken to be deZia-correlated in time corresponding to 
Markovian white noise 



{F,{t)F,{t')) =2D,, 5{t-t'). 



(7) 



The Dij are the diffusion coefficients in analogy with clas- 
sical Langevin equations. The diffusion coefRcients asso- 
ciated with Eqs. lO are calculated using the generalized 
fluctuation-dissipation theorem 6] in Appendix A. 
The Langevin equation for the signal field operator is 



04 = -764 - iRX^azu + Fa, 



(8) 



Here we included an heuristic field damping rate 7 to 
allow for atom-field interaction inside a cavity. From this 
equation, we can construct the equation 



A = -7(A - C) - ig^- 



(9) 



The effects of a thermal heat bath that accounts for a 
non-vanishing Fa^ are neglected by assuming low tem- 
perature. In this paper we will also assume two photon 
Raman resonance thereby set A = 0. 

We next define macroscopic variables by summing over 
the operators for all the particles (atoms/molecules) in 
the medium 



M 



N, 



b{c 



(10) 



and the population sum and difference operators Nt — 
Nb + Nc and N = Nb — Nc- The equations of motion for 
the collective atomic variables and the field are then 

M = -TM + gNA + Fm 
N = -TbiN + Nt) ~2g[A*M + M*A]+ Fn 
Nt = 

A = -jiA-O + gM, (11) 

with the total number of active particles {Nt) being con- 
served. The associated diffusion coefficients are shown in 
Appendix ^ 



Solution of these equations is facilitated by transform- 
ing them into equivalent c- number equations IS ISI; 3- 
This is achieved by putting all the operators in normal 
order chosen to be A^ , M^, N, M, A which establishes an 
unique correspondence between the quantum and classi- 
cal equations. Since the equations are already in normal 
order the operators are simply replaced with their classi- 
cal counterparts A* , M* ,Af, M, A. 

The c-number noise functions satisfy {J-{t)) — and 
{!F{t)T{t')) = 2T>5{t — t') just like their operator coun- 
terparts. The c-number diffusion coefficients will how- 
ever acquire additional terms arising from normal order- 
ing; all the non-vanishing coefficients are listed in Ap- 
pendix ^ The corresponding properties of the noise in 
the frequency domain are summarized in Appendix ^ It 
is worth mentioning here that the expectations in the c- 
number representation are now in the Glauber-Sudarshan 
P-representation of a thermal distribution which happens 
to be a Gaussian distribution of zero mean. This means 
that the Gaussian moment theorem applies and the first 
and second moments suffice to determine the distribution 
completely. 



IV. EQUATIONS AND ASSUMPTIONS OF THE 
MODEL 



We consider small fluctuations about the mean val- 
ues for both atomic and field variables Ai{t) = M.o{t) + 
5M (t) , Af{t) = Afo (t) + SAf{t) and Ait) = Aa (t) + 5A{t) . 
It is convenient to define two real variables T{t) = 
M*{t)Aoit) + Mit)A*oit) and S{t) = |^(t)P, with lin- 
earized fluctuations 

SS{t) ^ Ao{t)dA*{t)+AUt)SA{t) 

6T{t) = Ao{t)SM*{t) + A*o{t)SM{t) 

Trit) - Ao{t)TM*{t) + Al{t)TM{t)- (12) 

The relation to the properties of the generated field 
that we are finally interested in are as follows: The 
strength of the signal field, given by the number of gen- 
erated photons, and the atomic noise in the photon num- 
ber, quantified by its variance, are determined by 



ni{t) 
6nl{t) 



{al{t)a{t)) 
ni{t) 



IA(t)-C|^ 



\A,{t)\' 



{5S\t)). 



(13) 



To reduce clutter we will often leave out the expectation 
brackets being obvious from the context. The variance is 
obviously linearized and neglects relatively small terms 
quadratic in the correlations. We will primarily consider 
regimes where the generated field is much weaker than 
the input fields (as we justify numerically in Sec. IVII B|) 
and therefore we will replace |^oP — ICP where appro- 
priate and consistent. 
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adiabatic elimination of the averages of atomic 
variables 



Since the fluctuations have vanishing average, we can 
write separate equations for the mean and the Hnear fluc- 
tuations. The equations for the averages can be had by 
simply leaving out the Langevin forces 



dtMo = -TMo + gAfoAQ 
dtJVo = -TbiMo + Nt) ~ 2g[A*oMo 



M*oAo 



(14) 



Atomic/molecular transitions typically occur on a faster 
time scale than the variations in a radiation flcld, which 
means that the medium will adiabatically follow the 
field for weak coupling. Even when using fast pulses, 
the coherence generating pulses can be taken to be of 
sufficiently long duration for this to be true. There- 
fore taking averages over intermediate time scales, the 
time-derivatives in the equations for the averages of the 
medium variables Mo and Ao may be neglected in what 
is called adiabatic elimination |^ and we obtain algebraic 
relations for the coarse-grained averages: 



Nt 



1 



Mo 



C Ao 



(15) 



where we set \Ao\^ ~ as discussed earlier. We define 
some parameters here that we will use often and which 
will serve to simplify our expressions considerably: 



B = 



r.r' 



c = 



Nrg^ 

r 



(16) 



The parameter B is dimensionless and C carries the di- 
mension of inverse time. As we corroborate numerically 
later ^ 1, so the upper level is never strongly 

populated, A/q — —Nt and we can simplify by setting 
1 + ~ 1 in the rest of the paper. Substituting the 

atomic mean values we get an uncoupled equation for the 
average of the flcld variable 



dtAo ^ -liAo - - CAo. 



(17) 



fluctuations 



In describing the fluctuations, we cannot make argu- 
ments for adiabatic elimination of the medium variables 
as we did for the averages; this is due to the presence 
of the rapidly varying noise functions. An attempt to 
make such an approximation can lead to unphysical di- 
vergences in the correlations; we will discuss this issue in 
some detail in Appendix |D] Therefore we solve coupled 
equations for the fluctuations for both the fleld and the 
medium. We write the relevant equations in terms of the 
variables T and S 



dt6T = -rST^ 
dtSN = -Tb5N 
dtSS = -7(55 + 



2g5Af\Ao\'^ + gAfoSS + J't 
+ 2CSS - 2g6T + Tjsf 
gST. (18) 



In writing these equations, we have treated a-s a con- 
stant. We will see that in all the cases we consider in 
this paper this assumption is valid, because either Ao 
will have an unvarying steady state value or it will sat- 
isfy ^0 — ^ J even as a function of position when propaga- 
tion in free space is inovlved (see next subsection). It is 
clear from the equations above that the fluctuations can 
be fully expressed in terms of the correlations of the vari- 
ables T and Af. Using the results in Appendix^ we find 
that the two variables are uncorrelated with each other 
but the strengths of their autocorrelations are specified 
by the diffusion coefficients 



2Vtt - 2{T-n)NTB\^\^ 
2Vf^^ - ATi,NtB\^\'' 



(19) 



where we used the expressions for the averages in 
Eq. (|15|l . In writing these expressions we used the weak 
field approximation Ao — C and the assumption of weak 
upper level occupation ^ 1, and therefore in effect 

we can treat these coefficients as constants. 



C. Free space description 

The equations above are appropriate for describing ex- 
periments conducted in a cavity of damping rate 7. How- 
ever if we wish to describe experiments in free space we 
have to allow for field propagation through the medium 
and effectively replace 7 — *■ vdz-An approach discussed 
by Drummond and Carter llCll | outlined in Appendix 
10 leads to the appropriate space-time equation for the 
field variable 



idt + vd,)A{z,t)^gM{z,t). 



(20) 



This equation follows from Eq. ()C9(I derived in the ap- 
pendix, on assuming that all the input fields propagate 
collinearly in the z-direction and their amplitudes ai(2,3) 
vary little over the interaction length L, i.e. the length of 
the region over which the fields interact with the atoms, 
and therefore may be treated as constants. 

Thus for experiments in free space, equations for the 
field variable in (|17|l and IjlSI) are modified to 



{dt + vd,)Aoiz,t) = -CAo{z,t) 
{dt + vd,^)6S{z,t) = gST{z,t). 



(21) 



The equations for the atomic variables are formally un- 
changed, but they now carry both position and time argu- 
ment (z, t)] their spatial dependence arises through their 
interdependence on the field variable. 



V. PULSED INPUTS 

We first consider the case where the input laser fields 
are extremely short and fast pulses. In particular we 
take the pulses to arrive in sequence as is the case in 
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FAST CARS Q; the fields ai and a2 are allowed to 
interact first with the atoms for a duration Atj, creating 
the coherence, then after a delay Ai^ the field occurs 
for time Aia. This means that the atomic variables can 
now only depend on the first two fields, so we replace 
ylo ^ ^ in Eqs. (flTjl and ifTt^ and also set 5S = Q in the 
equations for the fluctuations of atomic variables: 



dtSAf = -TbS^^-2gST + T^r. 



(22) 



The main premise for this scenario is that when a rapid 
sequence of fast femtosecond pulses are used, the field 
scatters off the generated coherence before it has time 
to decay significantly. Thus for time intervals such that 
Atd + At3 <C r^^, the generated coherence may be 
taken to be constant. For longer durations however we 
have to allow for decay of the coherence A^o in Eq. (|15|l . 

During the signal-generating cycle when is present 
the mean coherence is described by 



dtMa = -TMo-giAo-O-^a- 



(23) 



After the coherence generating pulses arc turned off, the 
value of generated coherence from Eq. ifT^ A^o(O) = 
g£,Afo /F can be taken as the initial value of the coherence 
in the above equation (if Atd — 0) and since initially 
there is no signal field -4o(0) = It is then easy to see 
that the ratio of the magnitude of the first term to that of 
the second term is ^ |^|/|y^o — which is large for weak 
generated field and short duration pulses. Therefore we 
can neglect the second term and take the driving coher- 
ence to simply decay exponentially from the moment the 
generating fields are turned off Mo{t) = A4o(0)e^^*. 
The equation for the average field then becomes 



dtAa = -7(^0 -0- C^e 



-Tt 



(24) 



The fluctuations will likewise decay since the diffusion 
coefficients depend on the averages Aq and M.q. The 
diffusion coefficients in Eq. H19|) will then have an expo- 
nentially decaying time dependence, T>j^j^ ~ e"'"'* and 
T>rr ~ e~'"*. We note however that there is no con- 
straint on the duration A^c of the coherence generating 
pulses, so they can be taken to be sufficiently long to cre- 
ate a steady state value of the coherence before they are 
turned off, thereby justifying adiabatic elimination. 



A. pules: interaction in a cavity 

The mean value of the field variable after the third 
pulse Aq^AI-^ is easily obtained by integrating Eq. H24|) 



A(At3)-C--Cex 



-rAt3 



(25) 



The phase of is given by that of ^, since Aaj^ is real. 
When the duration of the third pulse is short Ata ^ F~^, 



A^[AH)-i^-CiAH. 



(26) 



The equal time fluctuations of the field are given by 

(552(At3) ^g^ dt dt'e-^(2At3-t-t')(5'r(i)(5r(i')l27) 
Jo Jo 

In order to evaluate this we first take a Fourier transform 
of the Eqs. (|22|l . the relevant definitions and properties 
of the fluctuations and noise in the frequency domain are 
described in Appendix IbI 

{r-^ioj)6T{uj) = 2g\^fdAf{uj)+TT{io) 
(Tb-iuj)dM{Lo) = -~2gdT{Lo) +Tj^{lu) (28) 

which we then solve to get 



[2g\^\^T^f + (Fb - ilu)Tt] 



(F — iuj)(Ti, — iuj) 



(29) 



Assuming weak excitation we set |1 4 [Y-iiJ \[ Th 
1. Using this expression for 5T{lo) we can evaluate the 
integrals in Eq. H27|l . The details of the calculations are 
shown in Appendix^ The general result result derived 
there is not particularly illuminating, instead here we 
consider the relevant limiting cases. In the short pulse 
limit A^s ^ F^^, we found that a Taylor expansion leads 
to an exact cancellation of the terms linear in At^ and 
we get a quadratic dependence on the interaction time 



B\^\^2V^^ , 2Vtt 

(r + Tb) 



F 



(30) 



A conspicuous feature of this limit is that 7 is absent 
in the expressions for both the mean and the fluctua- 
tion. Thus for pulses of duration shorter than the cavity 
damping time, the presence of the cavity has no effect 
on the signal or its associated noise due to the medium. 
This short pulse limit implies the hierarchy of time scales 

7-i»F-i,F,-i» Ai3. 

In the opposite limit of long pulse duration At^ 3> 
7~^ ^ F~^,Fj^^ the fluctuation is approximately 



SS^Ats) 



g^Atle-^^'^ 



rFfa7 



Bier 



2Prr 
F27 



.(31) 



Here we set the depopulation rate and decoherence rates 
to be equal F ~ F^, to mask unnecessary details and 
highlight the main feature, which is that both the mean 
and the fluctuations essentially vanish towards the end 
of a long third pulse as the driving coherence disappears. 
The conclusion is that, regardless of the duration of the 
third pulse, significant signal field is generated only dur- 
ing times satisfying < F^^,F^^. 



B. pulses: interaction in free space 

The last subsection showed that increasing the dura- 
tion of the interaction between the third pulse and the 
atoms will make a difference only upto a point, since we 
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are limited by the decay time of the coherence. In free 
space, this means that we will not get a stronger sus- 
tained signal by simply increasing the interaction length 
L. Therefore we will confine ourselves to sample sizes 
L < wF^^, iir^^, and we will use Eq. (|21l) where, taking 
the coherence to remain unchanged during the time of 
interaction, we set AQ{z,t) ^. 

A time- frequency Fourier transform and subsequent in- 
tegration over the interaction length gives 



{—ioj + vdz)Ao{z, ut) 
^ Ao{L,uj) 



-C£,5{uj 



(32) 



V 



and an inverse transform yields the delta- function S(t — 
L /v) that simply tells us that time is a redundant param- 
eter and the field can be specified by its position alone 



Ao{L) 



V 



(33) 



Likewise for the variance we consider the equations for 
the fluctuations in the frequency domain. Since we set 
Ao{z,t) ~ ^, the equations are given simply by 
Substituting the expression ST{u!) from Eq. (|29l) into the 
equation for the field fluctuation 



— iuj + vdzSS{z, uj) = gST{z, uj), 



(34) 



and integrating over the interaction length L we get the 
spectral density of the noise 



27ru2 (r2+Cj2)(r2+^2) 



The autocorrelation at equal times obtained by doing 
a partial fraction decomposition and an inverse Fourier 
transform defines the variance 



5S{Lf 



2t2 r 



(r + Fb) F 



(35) 



We see that both the mean and the variance are identi- 
cal to what we found in the short pulse limit when the 
interaction took place inside a cavity, bearing in mind 
that here L/v defines the time of interaction. This reaf- 
firms the conclusion that for short pulses there is no real 
advantage in using a cavity. 



VI. CONTINUOUS WAVE (CW) INPUTS 



A. cw: interaction in a cavity 

The average is found by integrating Eq. (|17|) , exactly in 
the form it is written, with the initial condition .4o(0) = 



(36) 



Unlike the pulsed case there are no obvious time con- 
straints, except those arising from possible damage to 
the sample by prolonged exposure to the fields, so we 
can take the long time limit t oo and we get a steady 
state signal 



1 



ii + cy 



(37) 



We note that if 7 <C C the generated field essentially 
vanishes, which suggests that we need to have 7 > C. 
In that case we can further use the assumption of weak 
signal field to conclude that 



— < 1 

7 



7 



(38) 



The fluctuations are best determined by writing the 
three coupled equations in the frequency domain in ma- 
trix form: 



(r-iuj) -2g\Ao\^ -gAfo 
2g {Tb~iuj) -2C 
-g (7 - iuj) 



Inverting the coefficient matrix Mp gives the solution 
[(Fb - iuj)gTT{oj) + 2.g2|Ap.^A^(a;)] 



5T 






6N 






SS 








5S{lo) 
det Mf 



det Mf 

(F — iuj){Tb — «w)(7 — iLo). 



{39) 



In computing the determinant in the denominator we in- 
troduced some simplifications using the arguments lead- 
ing upto Eq. (|38|) and we neglected the frequency depen- 
dence in the term proportional to -B|.4oP [refer to the 
comment following Eq. (|29|l ] . Squaring this gives the 
power spectrum (Appendix and then a Fourier trans- 
form and the adiabatic assumption 7 <i; F,F;, yields the 
steady state fluctuations 



We now consider the cases where all three input fields 
occur continuously and simultaneously. In this case 
steady state values for the medium variables may be con- 
sidered, and their decay does not put limits on the inter- 
action time as it did previously when using sequential 
pulses. But the calculation of fluctuations is complicated 
by the fact that the equations for the medium variables 
and the field variable do not decouple as they did for 
sequential pulses. 



27 



FFh 



F2 



(40) 



While the adiabatic assumption was not necessary it 
gave a simpler expression, the more general expression 
is shown in Appendix |E] What sets this case apart from 
the rest is that both the signal field and the variance in 
the steady state depend on the cavity lifetime, 7^^ and 
that dependence is linear. 
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B. cw: interaction in free space 

In order to find the signal strength, we do a Fourier 
transform of Eq (|21|l 

{-iuj + vd2)AQ{z,uj) = -CAo{z,uj) (41) 

and integrate it over the interaction length L to get 



iLuL/v 



A(0,c.)e-^^/". (42) 



As in the analogous case in a cavity, the time parameter 
is seen to be redundant after an inverse Fourier transform 
and we get simply 



Ao{L) - 



-CL/v 



(43) 



In calculating the fluctuations we take the Fourier 
transform of the equations for the atomic variables as 
they appear in Eqs. (|18|l and the field fluctuation as it 
appears in (|21l) . keeping in mind that all the fluctuation 
and noise elements will thereafter carry an argument of 
(z, uj). The three coupled equations yield an equation for 
the field fluctuation 



g [2g\^\^T^f + (Fb - ^c^)■Fr] 
(F — ia;)(Fb — iLu) 



. 9 [we + in ^^)gA/'o] .... 
with g{uj) ~ — — — 44 

(I — IUJ)(1 h — iw) 

We integrate this with respect to the spatial coordinate 
z, and thereby we get the spectral density 



1 L i_e2fle[eM]i 



27r V -2Re[g{uj)] 
4g^|^|'22?AAAf + ff^(rg + c.^)2I?rr 

(F2+L^2)(r2+^2) 



We use -B|^p <C 1 to simplify 



2Re[g{uj)] 



(45) 



We note this quantity achieves the largest magnitude 
when UJ = 0. In the case of short interaction length 
CL/v ^ 1 both the signal and the variance resemble 
those we got in the case of pulsed inputs in free space 



(46) 



SS'^{L,uj) 



1 [Ag^\^\^2V^^ + g^Tl + u;^)2VrT] 



27r 



(F2+c^2)(r2^^2) 



In the opposite time limit of long interaction length 
CL/v » 1, we find that the expression for the average 
value of the field variable Aq tends to become increasingly 
smaller. Equation (|43|l tells us 

n4(i) = |^o(i)-eP = |e(l-e-^^/")|2, (47) 



so that in this limit the input fields vanish and give way 
to the signal field. Our model however does not allow us 
to obtain an accurate expression for the variance in this 
limit. This is because the equations (|18|1 were based on 
the assumption that Ao does not change significantly and 
the generated field is relatively weak, and this is no longer 
true when L becomes large. Yet based on the fact that 
average value decreases over the length L and hence 
the diffusion coefficients also decrease, we could expect 
that the noise contribution from the active medium will 
actually be lower farther along the interaction length. 



VII. DISCUSSION OF RESULTS, INFERENCES 
AND ASSUMPTIONS 



We will now discuss in some detail the physical impli- 
cations of the results we obtained in the previous sec- 
tions, and also elaborate a bit more on our assumptions 
and some associated subtleties that we touched on while 
deriving those results. We begin by considering realistic 
numerical values for our parameters and variables and 
thereby provide concrete justification for the approxima- 
tions we made. 



A. Numerical Estimates 

For the purpose of numerical estimation we will use the 
specific example of an anthrax spore for which the Raman 
active molecule is dipicolinic acid (DPA) which consti- 
tutes 17% of the weight, the rest being mainly water. The 
number density of DPA molecules in an anthrax spore is 
Nt ~ 4 X 10^^ molecules m~'^, and the dimensions of the 
spore itself is 1 x 2 x 1 ^jlvo? . We will take all optical fre- 
quencies to be in the visible range v ^ 2-kc/\ ~ 4 x 10^^. 

Since some of the properties of DPA are not easily 
available, for those properties we will use the values for 
benzene, an organic molecule that has a similar struc- 
ture. Thus for instance we use the spontaneous Ra- 
man differential cross-section for Benzene ^ = 32.5 x 
10~'^^m2/sr/molecule, and we take the linewidths of Ra- 
man transitions in Benzene to estimate Fb ~ lO^^Hz. 

We first determine the strength of the Raman coupling 
R[j^ we can do that in two ways, assuming in either case 
that only a few intermediate states contribute: first using 
Eq. ^ , and taking /i ~ eag with og being the Bohr radius 



~ 4 X lO-^C VVJ's; (48) 



and secondly using the Kramers-Heisenberg formula for 
the differential scattering cross-section applied to spon- 
taneous Raman scattering 

'5~0.35x lO-^C^rnVj^s. (49) 
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So we will take \R[^\ ~ IQ-^ C^mVJ^s. Next we de- 
termine the field density and photon density from the 
expression for the field intensity (power/area) 



1 



1 



/ = -Wpeolf a| = -v{hv) — . 



(50) 



For strong lasers we could take the typical intensities to 
be / ^ 10^^ Wm^^, in which case we get 



\£a? ~ 5 X IQi'^N^C-^, 



n 
V 



V 



5 X 10^^m~-^(51) 



Taking all the input fields to have similar intensities and 
assuming natural linewidth Ft, ^ 10^^ Hz and decoher- 
ence rate F - lO^^ - iQi^ Hz, 



\£a\ 



FF,, 



10" 



(52) 



Indeed for this choice of parameters we are justified in 
taking the weak-excitation limit ^ 1, and in fact 

this will be valid till the field intensities increase to about 
lO^s Wm^^; we note that this is the intensity level at 
which cascade breakdown of air occurs at STP ll], so 
our weak- field assumption covers most realistic regimes. 
Thus we are quite justified in setting 1 -I- B\^\'^ ~ 1. 

Finally we address the issue of the interaction times, 
since many of our results assume short interaction times. 
In the case of pulsed inputs, F~^ sets the limits on the 
pulse duration to be Aia < 10~^^ s, which is certainly 
within the bounds of experimental capabilities using fem- 
tosecond pulses. In the case of propagation in free space, 
taking the dimension of an anthrax spore L ~ 10~^ m 
as the interaction length we find L/v ^ 10~^^ < F~^. In 
either case the constraints of short interaction times are 
likely to be satisfied for physical regimes of interest. 



B. Signal and noise 

The main observation we have from our calculations 
is that when the duration of interaction between the 
signal-generating field and the Raman-active medium 
At = A^s, L/v is short compared to F~^, the signal and 
the noise due to the medium have essentially the same 
theoretical description independent of the various exper- 
imental scenarios that we considered: 



714 



A^||i?43-Rl2p^^l?^2^^3A^2 



F2 

(r + Tfc) 



F 



n. X 



8(F - Fb) 
NrTb '■ 



(53) 



where we have used ^ 1 to set Ao — Nt and 

C ~ C; we also noted that the term involving Vj^/j^ is 
smaller by a factor of than the one involving Vtt 

and hence we only retained the latter. 



The exception to the above expressions was the case 
of continuous wave input in a cavity for which we got a 
steady state signal and noise given by: 



714 



5n. 



jV||i?43-Rl2p"-l»2?T-3 1 
F2 t2 



n. X 



87(r - Tb) 
NrVVb 



(54) 



The signal has the same form as the other cases with 
At 7^^. But the fluctuations differ by a factor of 
7/F. Since typically 7 <C F the noise noise will be less in 
this case, and because 7"^ is greater than the short time 
limits in the other cases the signal is bigger. In addition 
to this, the fact that there are no major constraints on 
the duration of the irradiations apart from their possible 
destruction of the sample, this scenario seems to be the 
best one from the signal to noise perspective. However 
that conclusion has to be moderated by the fact that 
when all three input fields happen simultaneously, there 
can be significant contributions from the non-resonant 
terms in our interaction Hamiltonian in Eq. particu- 
larly if we cannot achieve two-photon Raman resonance. 
This is where the sequential pulse scheme as in FAST 
CARS has an advantage. 

We can recast our expressions for the signal and fluc- 
tuations in terms of macroscopic variables. First using 
Eq. (O and the definition of classical field amplitude 
El — £iai we write the signal in terms of the classical 
field amplitudes: 



2eoer 



[X(3)]2|i?inij2|2|i?3|2x Ai2(55) 



where At — At^,L/v,^~^ depending on the experimen- 
tal configuration. Written this way we see that our ex- 
pression for the signal is consistent with what we would 
get if we integrated the classical equation Q for the 
slowly- varying amplitude with appropriate assumptions. 
A more practical representation would be in terms of the 
input intensities using Eq. (|50|l 



5ll 



It X 



c'^^^fhhhAt' 



8(F-Fb) 
TbNr 



for cw & cavity^ . (56) 



At this point we validate numerically the assump- 
tion of weak signal relative to the input fields. Using 
the parameters in the previous subsection we find the 
third order susceptibility to be x^^^ ^ lO"^'* C''N~^m~2. 
Taking all the input fields to have the same intensity 
li ~ Ii ^ I2 ^ I3 and the ratio of the signal intensity to 
the input field intensity is 



h 
h 



io-i°/Mi 



(57) 



For the short interaction times that we consider At ~ 
10~^^ s and input field intensities of li ~ 10^^ Wm"^ 
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this works out to be I4 ~ 10~^'^ li. This shows that even 
for higher input field intensities and increased density 
of active molecules, the signal field would still remain 
relatively weak. 



C. Comparison with Shot Noise 

Starting with the expression we derived for the gener- 
ated field, we can write the operator for the generated 
field in terms of the input fields: 



NT\R43Rl2\At ^ 

04 = = aial^as. 



(58) 



Then noting that 



a\a4 oc ni(n2 + l)n3 
a\a4a\a4 oc (nf + ni){nl + 3n2 + l)(n3 + n^) (59) 

and on assuming similar and large photon numbers in 
the input field, ni ^ n2 ~ ns = rij 1 we find that the 
variance corresponding to the shot noise is: 



x3xn^.3!i (60) 



Therefore the ratio of the medium noise and the shot 
noise is: 



(5n2 ^ 8 (r - Tb) rii/V 
BFb ^ Nt/V 



(61) 



Using the numerical values that we considered earlier we 
find 



10-2 - 10"^ 



(62) 



which shows that the noise due to the medium is less 
than the shot noise. The shot noise increases with in- 
creasing input field intensities while the medium noise 
increases with increasing number of particles. But in the 
ratio it is interesting that the behavior is exactly the 
opposite, the weight of the shot noise decreases with in- 
creasing intensity of the input fields while the weight of 
the medium noise decreases with increasing density of 
the medium. Therefore when the field intensity becomes 
very strong for a given medium density, shot noise could 
become lesser, but as we noted earlier the intensities can- 
not be much stronger than what we already considered 
without destroying the sample completely. On the other 
hand densities of active particles could be higher in other 
medium of interest. In the case of cw inputs in a cavity 
we have an advantage that the medium noise is further 
reduced by a factor of 7/r. 



VIII. CONCLUSION 

We developed a model based on Langevin equations 
which allowed us to get a purely analytic description 
of the signal and the quantum noise for Coherent Anti- 
Stokes Raman Spectroscopy. If we use sequential pulses 
(as in FAST CARS) or we arc close to two-photon res- 
onance the non-resonant terms would not be important 
and the quantum noise arising from the finite lifetimes 
and coherence times of the atoms/molecules would be 
a dominant source of noise. When interaction time be- 
tween the input fields and the medium is short, we found 
that the signal and the medium noise have the same be- 
havior in free space or in a cavity, and with pulsed inputs 
or with continuous waves. 

In particular we showed that if the driving fields do 
not vary much and the signal field is weak in compar- 
ison, the shot noise which represents the fundamental 
limits of noise is larger than the quantum noise due to 
the medium, and so the latter is not a limiting factor. 
Using a cavity to achieve steady state with continuous 
wave inputs leads to enhanced signal and lesser medium 
noise; this would be important close to two-photon res- 
onance and the non-resonant terms do not contribute to 
the background noise. 

Our calculations should be particularly relevant for 
novel experiments with newly developed femtosecond 
lasers and for fast spectroscopic characterizations of mi- 
croscopic agents in the air which could be organic ones 
like anthrax spores or inorganic suspensions or trace con- 
taminants. 

AlthoTigh our calculations in this paper were specific 
to CARS, the model we developed should be applicable 
to most coherent Raman schemes with minor changes. 
We did not resort to the commonly used adiabatic elim- 
ination when calculating the noise and we point out the 
significant errors that would arise from such an approxi- 
mation. Thereby we set the grounds for a more accurate 
understanding of quantum noise in stimulated Raman 
spectroscopy. In the regime of short interaction times 
we achieved a completely analytical description. But the 
Langevin equations we set up in our model describe a 
much broader physical regime, the solutions in general 
will have to be found numerically. 
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APPENDIX A: DIFFUSION COEFFICIENTS 

A quantum Langevin equations for an operator x(t) 
has the general structure 



(Al) 



with a deterministic part (not to be confused with the 
field operator defined in Eq. ©) and a stochastic part 
Fx- The diffusion coefficients 2Dxy — (F^Fy) associated 
with these equations are calculated using the generalized 
fluctuation-dissipation theorem |^ often called the 
Einstein relations: 



dt{xy) ^ (F^Fy) + (xAy) + (A^y) 



(A2) 



The non-zero diffusion coefhcients, corresponding to 
Eqs. © for the microscopic atomic variables, are: 

2D,,, = {2T^n){at) 

2D,,, = 2T{ac)+rb{ab), 

2D,b = -2D,c = Tb{a) 

2D,, = 2Dbb ^ -2Dbc = Tb{ab). (A3) 

From these using the definitions in Eq. (|10|l the non- 
vanishing diffusion coefficients for the macroscopic col- 
lective atomic operators M and N — Nb — are imme- 
diately obtained: 

2Dnn = ^TbiNb), 

2D^^,^ = {2T-Vb){Nb), 

2Dmn = 2Tb{M), 

2Dmm^ = 2T{N,)+Tb{Nb)- (A4) 

In transforming to c-numbers, by normal ordering, all the 
moments for the noise must remain unaltered. Gaussian 
noise is determined by the first two moments; the first 
moment (i.e. the mean) vanishes, the second moments 
must have the same time-evolution for a pair of operators 
i,y in normal-order and their c-number equivalents: 

dt{xy) = dt{xy). (A5) 

The fluctuation-dissipation theorem Eq. (|A2|) and its 
classical counterpart thereby relates the c-number diffu- 
sion coefficients to the operator diffusion coefficients: 

{T^Ty) = {F,Fy) + (xAy) + {A^jj) ~ (xAy) ~ ( Aj/) -(AG) 

Normal ordering of (xAy) + {Axy) gives some additional 
terms, and if x, y are not in normal order, we need to 
replace dt{xy) dt{yx) = dt{xy) — dt[x,y]. This de- 
termines all the non-vanishing c-number diffusion coeffi- 
cients for the collective atomic variables: 

2V^f^ = 2rb{NT+Af) -'ig{{M*A + A*M)), 
2Vm'M - [T^\T,){Nt+N), 
2Vmm = 2g{MA). (A7) 



APPENDIX B: NOISE SPECTRUM 

The Fourier transforms of the noise functions are de- 
fined by: 

F{uj) J dte"^'F{t)- F{t) = J duje'^^'F {ij).{Bl) 

For stationary processes the second moments or two- 
time correlations depend only on the time difference: 
{Fi{t)Fj{t -\- t)) — T{t). In the frequency domain this 
defines the spectral density P(w) of the noise: 



{F,{iu)Fj{uj')) = S{u; + oj')P,j{u) 



(B2) 



The two-time correlation and the spectral densities are 
time-frequency Fourier transforms of each other: 

P^M = ^ j dTe'^-Vir)- V{t)= j dre-'^^ P,j{lo). 

The fiuctuations considered in this paper are charac- 
terized either by (5-correlations or by exponentially de- 
caying correlations as a function of time difference for 
which the spectral densities are respectively constant and 
Lorentzian: 



r(T) = 2D,,5(t) 



r(r) 



-7|r| 



Pi. 



2A,— 

27 1 
72 ^2 27r ■ 



(B3) 



APPENDIX C: PROPAGATING FIELDS 

In order to describe fields propagating through a sam- 
ple of atoms in free space, we use the technique of Drum- 
mond and Carter 10]. Wc illustrate this technique by 
considering the propagation of the signal field V/^^ with 
wavenumber k in the direction of propagation (the z- 
axis); the part of the total Hamiltonian involving this 
field is 



H = hva\ai -t- Tig[dia'^ + H.a] 



(CI) 



For simplicity we have assumed a real interaction 
strength g =^ Ri^at^ = |i?43Q!3|. The interaction length 
L is divided into 2m + 1 equal segments of mean posi- 
tions zi — lL/{2m + 1) for / = — m, • • • ,m. Description 
of propagation requires multiple modes; a natural basis 
is provided by the normal modes for periodic boundary 
conditions on the quantization length L 



2'Kn 

IT 



(C2) 



with corresponding creation and annihilation operators 
cjj and Cn- Thus we can write the Hamiltonian as a linear 
combination of the sub Hamiltonians for each discrete 
segment summed over all the normal modes 



H.a. 



(C3) 



i.n.l 
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The index i accounts for the number of atoms in each 
segment, — A^T/(2m + 1). Introduce local operators 
for the slowly varying field amplitude in each segment 
through a discrete Fourier sum of the modes 



-(0 
04 



1 



n— — m 



J2m+ 1 

and the Hamiltonian can be written as 



(C4) 



m + la[ & 



t(0^t(0 



(C5) 



E™ 27rnc f 2TTin{l - 



n— — m 



(2m+l)L V 2m + 1 



The equation of motion for the slowly varying field am- 
plitude is then 



04 



(0 



(C7) 



We convert to c-numbers and take the limit of to — > cx3 
so that we have the following correspondence 



IL 



zi 



wiii^2m + 1 ^ c 
-I hm (2m + l) VaW)|,,^, 

m, — >r>Ci ' * 



2to + 1 
a/2to + la4 ^ 04(2, i) 

(C8) 



and thereby we arrive at the space-time dependent equa- 
tion of motion for the slowly varying amplitude of the 
propagating field 

^^ + ^-TAai{z,t)=gM{z,t). (C9) 

In the same way the space-time dependent noise opera- 
tors are defined by 

= lim (2TO+l)y /'^'(t) (CIO) 



and analogously their c-number counterparts. The space- 
time dependent c-number diffusion correlations 

{Tx{z,t)Ty{z,i)) 
= lim (2TO+l)2y(.F;'(i)^f 



APPENDIX D: PROBLEMS IN USING 
ADIABATIC ELIMINATION IN COMPUTING 
NOISE 

We pointed out earlier that adiabatic elimination of 
atomic variables cannot be applied to the fluctuations, 
we will now briefly discuss the consequences of doing so. 
We set 1 -I- B\^^ ~ 1 as justified earlier. We eliminate 
the atomic variables and write the consequent equation 
for the field variable inclusive of both the mean and the 
fluctuation 



A 



(Dl) 



Ta = - 



y-jV jT-^T 



In the case of a cavity, using the 8(t — t') correlation of 
the noise, we find 



J52(Ai3)=5- 



1 



-27At3 



27 



2I?rr , B\i\HVj^M' 



r2 



r^r 



While this reproduces the long term steady-state behav- 
ior correctly, the short-time behavior is linear in time, 
contrary to the quadratic dependence that we found ear- 
lier. 

The problem is more serious when we consider the free- 
space problem. We can understand it by considering the 
equation for the fluctuation in the field variable 



(D2) 



On integrating in the frequency domain and using the 
boundary condition M(0, i) = we find that 



\ n 
SA(L,t) = - / dzTA(z,t (L 



z)). (D3) 



As shown in Appendix|21the second moments of the noise 
are delta-correlated in space as well as in time 

{Ta{z. t)TA{z\ <')) = 2VAALf>{z ~ z')S{t ~ t') 

so that the second moments of the fluctuations of the 
signal at the output are 



{dAiL,t)5A{L,t')) 



2VAA^5{t-t'). (D4) 



At equal times this diverges, and therefore so docs the 
variance ^5*^, and this underscores why the adiabatic ap- 
proximation is inappropriate for describing the fluctua- 
tions. 



= L{T^{t)Ty{t))5{z~z'), (Cll) 
in the last step we used the correspondence 
(2m + 1) 



lim 



L 



-5w=5{z-z'). 



(C12) 



APPENDIX E: FLUCTUATIONS FOR 
INTERACTION IN A CAVITY 

We show the details of the derivation when the inter- 
action takes place in a cavity. From eqs. H27|) and H29|) 
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we get 
1 

^2^ 



At pAt 

g- I dt dt'e-^(2At-2t+r) 
/o Jo 



(El) 



with T = t—t'. On doing a partial fraction decomposition 
and carrying out the integrations 



(r2 



-27At 



-(rb+7)At 



r2) V27(7 + rb)r6 

1 _ „-27At 



g-27At _ g-(r+7)Ai 



+2X>rr 



27(7 + r)r (r2-72)r6 



V27(7 + r)r 



(r2 - 72)rb 



.(E2) 



In the hmit of short interaction time 7 At ^ 1 , it is easy 
to see that on doing a Taylor expansion of the exponen- 



tials as a function of 7At the zeroeth and first order terms 
cancel out. For example consider the coefficient of "IVrr 
in the last line which expanded to second order gives 



47^ At^ - (Ffc + 7)^At^ 

27(7 + rfc)rb " 2(rb - 7)(7 + rb)rb 



-272 Ai2 



Ae 



(E3) 



Thereby in this limit of short pulses we arrive at the 
expression for 88"^ {At) in Eq. (jSOJ. 

In the opposite limit of long time interaction times 
7Ai ^ 1 we get 



bS\At)^^- 



2Prr 



7(7- 



(r2 



r)r 
1 



1 



7(7 + rb)rb 7(7 + r)r 



(E4) 



In the case of sequential pulsed input we allow for de- 
cay of the coherence and the excitation and we get the 
expression in Eq. In the case of continuous wave 

inputs if we take take 7 <C F, Fb we reproduce Eq. H40|l . 
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